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Abstract 

We analyze the propagation of waves in unbounded photonic crystals. Waves 
are described by a Helmholtz equation with x-dependent coefficients, the scat¬ 
tering problem must be completed with a radiation condition at infinity. We 
develop an outgoing wave condition with the help of a Bloch wave expansion. 

Our radiation condition admits a uniqueness result, formulated in terms of the 
Bloch measure of solutions. We use the new radiation condition to analyze the 
transmission problem where, at fixed frequency, a wave hits the interface between 
free space and a photonic crystal. We show that the vertical wave number of the 
incident wave is a conserved quantity. Together with the frequency condition for 
the transmitted wave, this condition leads (for appropriate photonic crystals) to 
the effect of negative refraction at the interface. 

Keywords: Helmholtz equation, radiation, waveguide, Bloch analysis, outgoing 
wave condition, photonic crystal, transmission problem, negative refraction 
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1 Introduction 

Photonic crystals are optical devices that allow to mold the propagation properties of 
light. They usually have a periodic structure and are operated with light at a fixed 
frequency u. Due to their spectral properties (band gap structure), light of certain 
frequencies can travel in the photonic crystal, but, at other frequencies, the crystal 
is opaque. A large body of literature is available on this aspect of photonic crystals. 
Most contributions study a periodic medium, possibly with a compactly supported 
perturbation of the periodic structure. In contrast, we are interested in the interface 
between a photonic crystal and free space. 
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An interesting effect of such an interface is negative refraction. A recent discussion in 
the physical literature concerns the following question: Is negative refraction always the 
result of a negative index of the photonic crystal, or can negative refraction also occur at 
the interface between air and a photonic crystal with positive index? Our mathematical 
results conhrm the latter: The conservation of the transversal wave number can lead 
to negative refraction between two materials with positive index, as suggested in m 
In mathematical terms, the light intensity is determined by the Helmholtz equation 

— V ■ {a{x)'Vu{x)) = u{x), (1.1) 


which must be solved for m in a domain Q, u = u{x), x = {xi,X 2 ) € O. Here, we 
restrict our analysis to an unbounded rectangle O := M x (0, h) C but note that our 
methods can also be used in higher dimension, e.g. for H := M x (0, h 2 ) x (0, h^) C 


In (1.1), cj > 0 is a prescribed frequency and a = a(x) is the inverse permittivity of the 
medium. In an xs-independent geometry and with polarized light, the time-harmonic 


Maxwell’s equations reduce to (1.1) and u = u(x) is the out-of-plane component of the 
magnetic held. 

The coefficient a = a{x) describes the medium. We assume that the right half space 
{x = (xi,X 2 ) G Il|xi > 0} is occupied by a periodic photonic crystal with periodicity 
length £ > 0. Using the unit cube Y = (0,1)^ and the scaled cube = eY = (0,e)^, 
we therefore assume that the coefficient a = Y is W-periodic for xi > 0. We make the 
assumption that an integer number K of cells hts vertically in the domain, i.e. K = 
h/e G N. On the left half space {x = (xi,X 2 ) G H|xi < 0}, we set a = Y = 1. With 


a = Y and u given, problem (1.1) is an equation for u, but it must be accompanied by 
boundary conditions. 

We impose periodic boundary conditions in the vertical direction, i.e. we identify the 
lower boundary {x = (xi,X 2 )|x 2 = 0} with the upper boundary {x = (xi,X 2 )|x 2 = h}. 
In order to analyze scattering properties of the interface, we assume that the interface 
is lit by a planar wave. We consider, for a hxed wave-vector A: G the incident wave 


U.r 


_ 2'Ki'k-xle 


[X] = e 


( 1 . 2 ) 


To guarantee that Umc is a solution to (1.1) on the left, we assume = dvr^l/cp/e^. 
Since the Helmholtz-equation models a time-harmonic situation, we should think here 
of a solution of the wave equation in the form f/inc(a;,t) = Uinc(a^)e“*‘^* = exp(i[27rfc ■ 
x/e — ut]). We always consider fci > 0 such that U\-ac represents a right-going wave. In 
addition, we assume that the incident wave respects the periodicity condition in vertical 
direction, i.e. \ equivalently, k 2 h G eZ. 

With the incident wave Ujnc at hand we can now describe - at least formally - the 


boundary conditions for solutions u of (1.1) as Xi — )■ ±cxd. We seek for u such that (i) 
u satishes an outgoing wave condition as xi —)■ oo and (ii) u — Uinc satishes an outgoing 
wave condition for xi —>■ — oo. This leads us to our first question: 


Question 1: How can we prescribe radiation conditions in periodic media? 

The answer to Question 1 is intricate and requires a detailed study. We will use Bloch 
expansions and Bloch projections to formulate our new outgoing wave condition in 
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Figure 1: The geometry of the transmission problem for iF = 10 (number of cells in 
vertical direction). An incoming wave hits the boundary of a photonic crystal. We are 
interested in the waves that are generated in the photonic crystal. 


Dehnition |3.3 
subsection. 


In order to motivate our choice, we sketch some background in the next 


Once we have a precise formulation of the scattering problem, we can turn to the 
application: What can be said about the transmission problem? When an incident wave 
f/inc lights the interface, it creates waves inside the photonic crystal. These waves are 
described by u on {xi > 0}, our aim is to characterize u. Performing a Bloch expansion, 
we write u as a superposition of Bloch waves. In this superposition, we expect that there 
appear only waves that satisfy two requirements: (a) the Bloch frequency coincides with 
the frequency uj. (b) the vertical wave number of the Bloch wave is k 2 (“conservation 
of the vertical wave number”). 


Question 2: Let u be the solution of the transmission problem for the 
incoming wave 17inc- Does the Bloch expansion of u on the right respect the 
frequency condition and the conservation of the vertical wave number? 

A positive answer to Question 2 provides information on the negative refraction phe¬ 
nomenon. The requirements (a) and (b) are used in [27] to explain negative refraction 
without referring to a negative index material: Denoting the m-th Bloch eigenvalue for 
the wave-vector j E Z := [0,1]^ as the photonic crystal can have the property 

that the three conditions (a) fio{j) = (b) j 2 = ^ 2 , and the additional condition (c) 

ei ■ Vj/io(j) > 0 (the group velocity has a positive Xi-component), determine j uniquely. 
For an appropriately chosen held a, an appropriate frequency u and an appropriate in¬ 
coming wave vector k, we have the following situation: 62 ■ k = k 2 is negative, but 
the solution j satishes 62 ■ Vjp^oij) > 0. This means that a light beam that hits the 
interface from above (/c 2 negative in free space implies that the vertical group velocity 
is negative) produces a light beam in the photonic crystal that is directed towards the 
top (vertical group velocity is positive, 62 • Vjfio{j) > 0). With this mechanism, the 
conditions (a)-(b) can lead to negative refraction. This is outlined in [27|, where a 
specihc photonic crystal is described and the negative refraction effect is supported by 
numerical results. We note that a quite different interpretation is given in na. 
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We will answer the above Questions 1 and 2. The precise answers are more complex 
than one might expect at hrst sight (we sketch some of the principal difficulties in the 
next two subsections). We show that our outgoing wave condition of Definition 3.3 is 


reasonable by proving a uniqueness result: Theorem 1.3 yields, in a weak sense, the 
uniqueness of solutions in terms of the Bloch measure. Question 2 is answered with 
Theorem ]_A: If m is a solution that satishes the outgoing wave condition, then the 
corresponding Bloch measure is concentrated in those frequencies that respect (a)-(b). 


The mathematical description of our results is given in Section 1.3 


1.1 Outgoing wave conditions 

Although we use our results to analyze negative refraction, the core of our mathematical 
theory is more general: We develop an outgoing wave condition for the Helmholtz 
equation in a periodic medium. In this section, we sketch some background concerning 
radiation conditions, mainly in free space. Our aim is to demonstrate the importance 
of radiation conditions, to show the intimate link between radiation conditions and 
uniqueness results, and to motivate our mathematical approach. 


The Helmholtz equation (1.1) has been studied already by Euler and Lagrange, 


but Helmholtz was the hrst who expressed solutions in bounded domains with a rep¬ 
resentation formula [IH]. In unbounded domains, one faces the problem of boundary 
conditions at inhnity. We recall that two fundamental solutions of the Helmholtz equa¬ 
tion for X G are given by 


^out(^) • 



and 

ffinc(a:) := 

X 


X 


(1.3) 


With the time-dependence the solution Uout represents an outgoing wave. Mine an 

incoming wave. Outgoing waves are expected to be the building stones of solutions of 
scattering problems, incoming waves should not be present in the expansion of solutions. 

Sommerfeld introduced in EH for dimension n = 3 a radiation condition; until 
today, it is the standard outgoing wave condition in free space and is named after him: 


\x 


— iuju){x) —)• 0 as |x| 


—>■ oo 


(1.4) 


The solution Wout satishes (1.4) and is therefore admissible. Mine does not satisfy (1.4) 


and is not accepted as a solution. Sommerfeld justihed his radiation condition with 
a uniqueness proof: Prescribing boundary data on an obstacle (the scatterer) and the 


radiation condition (1.4) at inhnity, the Helmholtz equation has at most one solution. 


Actually, Sommerfeld demanded two further properties to guarantee uniqueness, but 
the results of Rellich (today known as “Rellich Lemma”) showed that the additional 
assumptions are not necessary [Mj, see [35] for the historical background. 

For two reasons, we cannot use the Sommerfeld radiation condition. The hrst is 
that we consider x-dependent coefficients a. The interest in x-dependent coefficients 
is not new: Sommerfeld himself studied the case that a takes two diherent values in 
two disjoint half-spaces, Jager studied in [21] coefficients a that stabilize to constant 
coefficients for |x| —)■ oo. Our situation is diherent, since a is periodic in the right 
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half plane. The second reason is that we study a waveguide such that, in the above 
sense, our situation is neither one- nor two-dimensional. For constant coefficients, the 
elementary solution in a strip M x (0,1) is with k = (ki, K 2 ) and \k\^ = u, which is 
right-going for ki > 0 and left-going for ki < 0. The solution has no decay (similar to 
the one-dimensional case), but the expression dx^u — ioju does not vanish for right-going 
waves due to the presence of k, 2 - 

The idea of our outgoing wave condition is simple: Every function on a rectangular 
domain can be expanded in Bloch waves. We demand that the expansion of the solution 
contains only outgoing waves. The precise form of the outgoing wave condition is slightly 
more technical since we have to consider restrictions of the solution u to large rectangles 
(in order to have a small contribution of non-periodicity effects). 

We emphasize that, even though the conditions become more technical, we follow 
the historical pathway: The expression in (1.4) can be understood as a projection of 
the solution u onto incoming waves. Outgoing waves are hltered out and it is demanded 
that the remainder is small for large radii. Our outgoing wave condition (3.6) is: At 
the far right, the solution u can be expanded in a Bloch series that contains only right¬ 
going waves. With this requirement, we follow once more Sommerfeld who writes in 
|37j : “at inhnity u must be representable as a sum (or integral) of waves of the divergent 
traveling type.” For an extensive study of homogeneous media we refer to [Hi- 


On radiation in waveguides and photonic crystals 

An important contribution is the recent work of Fliss and Joly [16], which is also con¬ 
cerned with outgoing wave conditions and the existence and the uniqueness of solutions 
for periodic wave-guides. Essentially, the outgoing wave condition of [IB] (for xi —)■ -|-cxd) 
reads 

u{.x)= ^ alU^{x) +w^{x), (1.5) 

where Niu) is a hnite index set, are real coefficients, are right-going Bloch waves 
and w~^{x) is exponentially small for xi —)■ -|-cx). The setting of the problem differs in 
one important point from ours: [I6] studies a medium which is identical at the far left 
and at the far right, which allows to use global Floquet-Bloch transformations; this is 
not possible in our setting. Below we give a more detailed comparison of our results to 
those of [TB] . 

Another radiation condition in a waveguide with varying index in transversal direc¬ 
tion is studied in [6]. The “modal radiation condition”, formulated in Dehnition 2.4 of 
[6], demands for solutions u of the Helmholtz equation that 

{fFu{x, .))(A) = for ±x > a 

holds for every A. Here, fF denotes a generalized Fourier transform and x is the longi¬ 
tudinal independent variable. As in our radiation condition, it is demanded that only 
outgoing waves instead of e+^l*l) are present. The new feature of our approach 

is that it covers media with oscillations also in longitudinal direction. We cannot use 
methods that rely on separation of variables and the Fourier transform must be replaced 
by a Bloch transform. 
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Also in |1], the radiation of waves inside a photonic crystal is investigated, and 
the setting also uses an interface between a photonic crystal and free space. The 
fundamental difference to our work is that in [1] the underlying frequency uj is assumed 
to he in a band-gap of the photonic crystal. For this reason, waves in |3] are found to 
decay exponentially in the photonic crystal and no explicit radiation condition must 
be formulated. We mention [30] and the references therein for other approaches to 
radiation conditions, also based on Poynting vectors and incoming and outgoing waves. 

For numerical calculations, one is interested in replacing the unbounded domain by 
a bounded domain. In this case, one asks for appropriate boundary conditions that 
must be imposed on the boundary of the bounded domain. This point of view leads 
to the construction of Dirichlet-to-Neumann maps or similar ideas [HttSlE]. Other 
key-words are perfectly matched layers [23| or transparent boundary conditions. These 
approaches may be an alternative to the outgoing wave condition that we suggest here. 
However, we are not aware of any result for such boundary conditions that implies our 
uniqueness theorem. 

1.2 Uniqueness and negative refraction 

Following Sommerfeld’s example, we accompany our outgoing wave condition in pho¬ 
tonic crystals with a uniqueness statement. With this result, we can treat the applica¬ 
tion on negative refraction. Several problems must be tackled in this process and the 
uniqueness result is, unfortunately, not as strong and simple as one would like it to be. 

It is an essential feature of the Helmholtz equation that, even without source terms 
and with homogeneous boundary conditions, solutions may be nontrivial. One example 
is the bounded domain O = (0,1) C with the solution u{x) = sin(7ra;) for cu = vr. A 
more relevant example in higher dimension (2 or 3) is the Helmholtz resonator: When 
UJ coincides with the resonance frequency, there is a nontrivial solution to homogeneous 
boundary conditions, see [36]. For regular exterior domains, the Sommerfeld condition 
implies uniqueness: The Helmholtz operator has only a continuous spectrum and no 
point spectrum. We emphasize that this is true only for the Helmholtz equation with 
constant coefficients. 

In our case of non-constant and (looking globally) non-periodic coefficients, there 
can be nontrivial solutions to the homogeneous Helmholtz equation (satisfying also 
a radiation condition). In general, such solutions can be localized modes or waves 
that travel vertically. An example for the hrst (corresponding to a point spectrum of 
the operator) are standing waves in a photonic crystals with a point defect, compare 
e.g. [22], Chapter 5. In the case of a line defect (or in our situation of an interface 
between free space and photonic crystal), one expects nontrivial solutions travelling 
along the interface, see e.g. [281 1^ . Concerning the mathematical analysis of defects 
in a photonic crystal and the possibility that they support modes (and hence act as a 
waveguide) see [51 US]- In a vertically periodic setting it was shown in [20] that a line 
defect cannot support bounded modes. 

The strong uniqueness result of [TB] is htting in this background: In a situation where 
the surrounding medium is perfectly periodic, the radiation condition of [16] implies 
uniqueness for non-singular frequencies. Vertical waves in the crystal are excluded 


Outgoing wave conditions and transmission at interfaces of photonic crystals 


7 


by the non-singularity assumption, waves along a line defect and localized waves are 
excluded by the absence of defects. 


Instead, our uniqueness result for the transmission problem must deal with the fact 
that the interface can support nontrivial solutions. Furthermore, we want to admit also 
singular frequencies (allowing for vertical waves in the crystal). Our uniqueness result 
states: Imposing the new radiation condition in photonic crystals, for non-singular fre¬ 
quencies, every homogeneous solution has a vanishing Bloch measure. Loosely speaking: 
the solution vanishes far away from the interface. For general frequencies, the radiating 
solution may contain vertical waves. See Theorem 1.3 for both results. 


A more technical problem will accompany us along the way to a radiation condition 
and to the uniqueness result: The geometry is not globally periodic and the solution u is 
not periodic (and u is, in general, not periodic on any rectangle in the right half plane). 
For this reason, neither a Bloch transformation of u nor a periodic Bloch expansion of 
u are meaningful as such. We will have to truncate u on large squares at the far right 
and consider the Bloch expansion of the result. We must use large squares in order to 
achieve that the truncation process introduces only small errors. 

Bloch measures (as used e.g. in [1], pp. 182-183) are the appropriate tool for the limit 
analysis, which is necessary for the following reason: A periodic Bloch expansion uses 
a discrete set of frequences j. In general, not even the elementary frequency condition 
ho{j) = (fhe Bloch wave frequency coincides with the frequency of the Helmholtz 
equation) can be satisfied in a discrete set of frequencies j. For this reason, we cannot 
expect that the Bloch expansion of u (at a finite distance) satisfies conditions such as 
po(j) = Instead, we must introduce a limiting object (the Bloch measure). Our 
aim is to derive properties of this limiting object. 

Regarding other mathematical approaches to related problems, we mention [2113], 
where the authors investigate diffraction effects in time-dependent equations. In [1], 
the spectrum of an elliptic operator in a periodic medium is investigated. We use some 
methods of [T] , in particular in the pre-Bloch expansion. Moreover, the above mentioned 
problem of waves that are concentrated at the interface of the photonic crystal has a 
counterpart in [1]: The part of the spectrum that is related to the boundary layer 
cannot be characterized explicitly (in the sense of [1], where the sequence £* —)■ 0 is 
fixed, and in contrast to |9|, where the sequence —)■ 0 is chosen appropriately). 


We close this section with more references to negative refraction effects. Negative 
refraction can be a consequence of a negative index material, see [31] for the effect and 
[3 El sni Ea ES] for rigorous results, obtained with the tools of homogenization theory. 
In [T^ [32] . the negative refraction effect is explained in the spirit of negative index 
materials. But negative refraction can also occur without a negative index material, 
see [EZ]- We note that the photonic crystals in [TE] and in [E7] are identical and that 
they do not have a negative effective index in the sense of homogenization. With the 
work at hand we support the line of argument of |27|. 
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1.3 Main results 


Throughout this article we consider the following parameters as hxed: The frequency 
(u > 0, the height h > 0 of the waveguide, the periodicity length £ > 0 with K = 
h/e eN, and the wave number /c G of the incident wave with k 2 h G eZ, ki > 0 and 
The underlying domain is hi := M x (0, h) and the coefficient held is 
a = : n —)■ M. We assume 0 < a* < a{x) < a* < oo Wx E Q and a = 1 on {xi < 0}, 

but the latter assumption is not essential. We demand a E with e-periodicity with 
respect to Xi and X 2 on {xi > 0}. 

We use Bloch expansions of the solution. Let us give a description of our results, 
where the superscript “±” indicates that we study ±xi > 0. The Bloch expansion uses 
two indices, m G No = {0,1,2,...} numbers the eigenfunctions in the periodicity cell 
and the Bloch number j E Z := [0,1]^ measures the phase shift along one periodicity 
cell. We collect the two indices in one index as A := (j, m) E I := Z x Nq. To every 
A G / we associate a Poynting number E M, see (3.1). For the Bloch wave with 


index A, the number is a measure for the hux of energy in positive Xi-direction. 
We introduce the outgoing wave condition (on the right) 


|nto(4)r^o 


as 


i? —)■ OO 


( 1 . 6 ) 


RY^ 


Here u~^ is, up to periodic extensions, the function u~^{xi,X 2 ) = u{Re + xi,X 2 ). The 
map n^Q is a projection onto those Bloch waves that correspond to an energy flux to 
the left (i.e. incoming waves, P^ < 0). The precise description is given in Dehnition 
3.31 The outgoing condition on the left is analogous; Bloch waves that correspond to 


an energy flux to the right (i.e. incoming waves, P^ > 0) are excluded on the far left. 

Our results are formulated with the help of index sets. Waves with vertical energy 
flux (or no energy flux) correspond to A G / := Z x No in 

/i:={Ae/|P^ = 0} , 

and, for a given m E Nq, to j E Z = [0,1]^ in the index set 


■= {j ^Z\ = o| = {j G Z I (j, m) Elto} ■ 

The statements below are meaningful for general frequencies a; > 0. Unfortunately, 
we are only able to prove theorems for moderate frequencies, as expressed in the follow¬ 
ing assumption. It demands that the frequency of the wave is below the energy band 
corresponding to the index m = 1 (below the second band). 

Assumption 1.1 (Smallness of the frequency). We assume on the coefficient a and 
the frequency oo that 

jGZ, m>l 

and /^m(j); where are the Bloch-eigenvalues. 

Our results concern solutions u of the transmission problem, specihed as follows. 




Outgoing wave conditions and transmission at interfaces of photonic crystals 


9 


Problem 1.2 (Transmission problem). We say that u G solves the scattering 


problem if it satisfies the Helmholtz equation (1.1) in Vt = M x fi)^h) with h = Ke 
and periodic boundary conditions in the X 2 -variable. We furthermore assume that it 
is generated by the incoming wave Uinc of (1.2) in the following sense: u satisfies the 
outgoing wave condition (3.6) on the right and the difference u — f/inc satisfies the 


outgoing wave condition (3.7) on the left 


Onr nniqneness resnlt characterizes the Bloch measnres of a difference of two 


solntions (the Bloch measnres are introdnced in Dehnition 4.2). The theorem below 


yields that, for large valnes of |a;i|, the difference of two solntions does not contain 
Bloch waves with an eigenvalne index larger than 0. Fnrthermore, only those waves 
can appear that satisfy all of the following three reqnirements: They correspond to 
the imposed freqnency u, they correspond to vertically periodic waves, they transport 
energy in vertical direction. 


Theorem 1.3 (Uniqneness). Let Assumption 1.1 on the frequency u be satisfied. For 
the incoming wave f/inc with wave vector k, let u and u be two solutions of the trans¬ 
mission ProblemlL^ For I E Nq, let be the Bloch measures that are generated by 

-u. Then: 


the difference v := u 


L.OO 


= 0 for / > 1, 


supp(!/J^) c {j e 21 pj(j) = e z/A'} n ^ 


( 1 . 8 ) 

(1.9) 


An immediate conseqnence of Theorem |1.3 is the following resnlt for freqnencies uj 
that do not snpport vertical waves. 

Corollary 1.4 (Uniqneness for non-singnlar freqnencies). Let the situation be 
as in Theorem M and let the frequency u be non-singular in the sense that 
{j G Z\ /ig (j) = , j 2 ^ K'^ f] J=o,o = 0- Then the difference v := u — u of two 
solutions of the transmission Problem L2 has a vanishing Bloch measure. 


Onr second main resnlt shows that the transmission of an incoming wave occnrs 
in snch a way that two qnantities are conserved: The vertical wave nnmber and the 
energy. 


Theorem 1.5 (Transmission conditions). Let Assumption 1.1 be satisfied, let k be the 
wave vector of the incoming wave f/inc- Let u be a solution of the transmission problem 
with I E 'Nq, be the Bloch measures that are generated by u. Then 


1.2 


and let 


moo ~ 0 ^ 1 


supp(i'too) C {i € Z I pj(j) = , ji € Z/A} n ({j € Z I it = kt] U Jf(, „) . 

( 1 . 10 ) 

As above, we have the following corollary for non-singnlar freqnencies. 

Corollary 1.6 (Transmission condition for non-singnlar freqnencies). Let the situation 
be as in Theorem \1.5 and let the frequency u be non-singular. Then the Bloch measure 
^Ooo satisfies 


supp(i^Joo) C [j E Z\ /i^(j) = and ja = ^2} • 


( 1 . 11 ) 
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1.4 Further comments on the main results 

On the uniqueness result. We recall that we expect the existence of solutions 
that are supported on the interface between photonic crystal and free space. For this 
reason, uniqueness results can only provide information “far away from the interface”, 
i.e. information on the Bloch measure. 

Our results are 


A weakness of our uniqueness results concerns Assumption 1.1 


proven under the assumption that the underlying frequency u is in the hrst band (more 
precisely: below the second band). Our conjecture is that our uniqueness result remains 
valid for arbitrary frequencies, stating that supp(r')^^) C {j \ ^f{j) = G Z/iF} O 

I for every I > 0. Due to a lack of orthogonality properties in the sesquilinear form 
b (see Section]^, we must exploit the frequency assumption in our uniqueness proof. 


Relations to Fliss and Joly jlffj/ . The contribution [11] contains strong results: 
1. A uniqueness result in the classical form (due to the absence of an interface that can 
support waves and due to the restriction to non-singular frequencies). 2. An existence 
result, based on a limiting absorption principle. We note that also the existence result 
of [16] uses global Floquet-Bloch transformations and is therefore not easily adaptable 
to our setting. We remark that our outgoing wave condition is weaker than the one of 
This means that, apart from the problems due to the non-periodic 


see Lemma 3. 


geometry, an existence proof should be simpler for our outgoing wave condition. 

We mention at this place that our outgoing wave condition differs in one point with 
all existing radiation conditions: Our condition does not use explicitely the frequency 
u. We regard this as an advantage: our condition might be applicable also in time- 
dependent problems. 


A possible scaling in e > 0. In all our theorems we keep the length scale e > 0 
hxed. In other words: the wave-length 1/uj and the periodicity length e are both of 
order 1. It is very interesting to analyze the behavior of light in small micro-structures, 
i.e. to analyze the limit e —?• 0. The limit can be performed in two settings: In the 
classical homogenization problem, one keeps u (and hence the wave-length) hxed and 
analyzes the behavior of solutions u = as e ^ 0. This approach was carried out 
e.g.in [ZlIHlIiniEslEe]. 

The second setting regards the limit £ —)■ 0 in a situation where the wave-length 
of the incoming wave is also of order £. This is the scaling that is suggested by our 
notation in (1.2), which corresponds to a frequency oj = oj^ = Loosely speaking, 

our Theorem 1.5 yields in this scaling: Solutions to the scattering problem with 
incoming wave (1.2) for hxed k consist, at a hxed distance xi > 0 from the interface 
and in the limit £ —)■ 0, only of Bloch waves that correspond to the frequency u:^ and 
to the wave number ^2 (up to vertical waves). 


Outline of this contribution. Bloch expansions are described in Section In 
Section we dehne energy hux numbers and corresponding index sets; these are used 
to dehne the new outgoing wave condition. In Section we dehne Bloch measures. 
Theorems 11.31 and 11.51 are shown in Section 14.31 
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2 Bloch expansions 

2.1 Pre-Bloch expansions 

We start our analysis with a discrete expansion. This discrete expansion is the first 
stage of a Bloch expansion and closely related to the Floquet-Bloch transform. We 
apply it to the h-periodic function u(xi, ■). The subsequent result appears as Lemma 
4.9 in [T]. 

Lemma 2.1 (Vertical pre-Bloch expansion). Let K G N be the number of periodicity 
cells and let h = eK be the height of the strip M x (0, h). Let u G x (0, h); C) 

be a function. Then u can be expanded uniquely in periodic functions with phase-shifts: 
With the finite index set Qk '■= ■ ■ ■, 

u{xi,X2) = ^ (' 2 . 1 ) 

j2&QK 

where each function <hj 2 (a;i,-) is e-periodic. The equality (2.1) holds in X 

(0,h);C). 

Sketch of proof. We sketch a proof (different from the one chosen in [T]), considering 
only u = u{x 2 ) and h = 1. Expanding u in a Fourier series, we may write 


u{x2) = ^ f3k2 


glTxik2X2le 


k2Ci£L 


For every j 2 £ with j 2 < 1 (he. for every j 2 ^ Qk) we set 


®i2(a^2) := 

k2&32 +1- 


fdkn e 


2m{k2-n)x2le 


With this choice, each is e-periodic and (2.1) is satished. 


( 2 . 2 ) 


(2.3) 


□ 


For the above pre-Bloch expansion we dehne the projection on a vertical wave num¬ 
ber /c 2 as follows. 

Definition 2.2 (Vertical pre-Bloch projection Ll^®'^*). Let u E x (0, h);C) with 

h = eK be a function on a strip and let ^2 G Qk be a vertical wave number. Then, 


expanding u as in (2.1), we set 

ng'‘«(ii,X2) := 


(2.4) 


The projection is an orthogonal projection: For e-periodic functions <F and <F an d 


A.l 


indices k 2 ^ ^2 there holds dx 2 = 0 by Lemma 

We will later use the following fact: If u is a solution of the scattering problem with 
incident vertical wave number k 2 , then also the projection is a solution of the 

scattering problem. Together with a uniqueness result for solutions, we can conclude 
from this fact that the vertical wave number is conserved in the photonic crystal. 

Below, we have to deal with the following situation: For a function u on a strip with 
height h, we can perform a pre-Bloch expansion. We may also extend u periodically in 
the vertical direction and perform a pre-Bloch expansion of the extended function on a 
wider strip. We hnd that both constructions yield the same result. 
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Remark 2.3 (Vertical pre-Bloch expansion of a periodically extended fnnction). Let 
K = h/e G N denote the number of periodicity cells in vertical direction and let u G 
X (0, h)) he a function with vertical pre-Bloch expansion 

u{xi,X2)= ^ 
j2&QK 

Let R E N be a multiple of K and let u be the periodic extension of u to the interval 
{0,eR) in X 2 -direction. Then u E L'LM X (0,e-R)) has the vertical pre-Bloch expansion 

u{xi,X2)= ^ (2.5) 

j2&QR 

where the coefficients according to the finer grid Qr satisfy 

(a:) = | ° ^fj2^QK, 

^ \ $.^(a;) ^fJ2 eQK. 

The statement follows immediately from the nniqneness of the pre-Bloch expansion. 
Remark |2.3| explains onr choice concerning scalings: Given a seqnence of fnnctions ur, 
dehned on a seqnence of increasing domains, at hrst sight, one might hnd it natnral 
to rescale ur to a standard domain and to analyze the seqnence of rescaled fnnctions. 
Instead, we work with the seqnence ur on increasing domains. In this way, one index 
j E Z always refers to the same elementary wave, which allows to investigate the Bloch 
measnre limit. 

Pre-Bloch expansion in two variables. For a fnnction u that is dehned on a rect¬ 
angle and that is periodic in both directions, the pre-Bloch expansion in two variables 
can be dehned by expanding hrst in one variable and then in the other. 

For fnnctions m on M x (0, h) the sitnation is more difhcnlt, since u is not periodic 
in xi-direction. In order to expand in both directions, we trnncate u with a cnt-oh 
fnnction p : M x [0, h] with compact snpport. For convenience, we assnme that the 
snpport of rj is contained in the square [0, h] x [0, h]. 

The truncation of u is dehned as w{x) := u{x)r]{x). We expand w (on the square 
[0, h] X [0, h]) in both directions in a pre-Bloch expansion, using the vector j = (ji, 72 ) ^ 
Qk X Qk and x = {xi,X 2 )-. 

w{x)= (2.6) 

J&QrxQk 

The functions now ^-periodic in both variables. Due to orthogonality 

there holds {h = eK) 

where -jy \^jf jy denotes the mean value. 
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2.2 Bloch expansion 

With the help of the pre-Bloch expansion we construct now the Bloch expansion. This 
step consists in developing each of the periodic functions for j = ( 71 ,^ 2 ) in terms of 
eigenfunctions of the operator 


:= - (V + 2Tiij/e) ■ (V + 27rij/e)) . 


(2.7) 


The operator acts on complex-valued functions on the cell W with periodic boundary 
conditions. It appears in the analysis of (1.1) for the following reason: Let be an 
eigenfunction of with eigenvalue then there holds 


We see that is a solution of the Helmholtz equation on the right half-plane 

if and only if = o;^. 

We have to distinguish between xi > 0 and xi < 0. On the right, the expansion 
is performed with as above, with the periodic coefficient a'^ = a^{x). On the left, 
expansions are performed according to = 1 with the operator Cj := — (V -|- 27rij/e) ■ 
(V -I- 27rij/e). The result is a classical Fourier expansion of the solution. 


Definition 2.4 (Bloch eigenfunctions). Let j G [0,1]^ be a fixed wave vector, 
denote by family of eigenfunctions of the operator of (2.7). 


We 

The 


labelling is such that the corresponding eigenvalues pfn{j) are ordered, fim+iij) > Tm{j) 
for all m E No- Similarly, (^ym)meNo family of eigenfunctions of the operator 

Cf and p-fiU) corresponding eigenvalues. We normalize with -jy = 1- 


A standard symmetry argument yields that, after an appropriate orthonormaliza¬ 
tion procedure for multiple eigenvalues, all functions T with j E Qk are 

orthonormal in the space L‘^{KYy<C) = Lf{KYy<C) (the sharp symbol is sometimes 
used to indicate that one thinks of periodic functions, but, of course, in the case of L^, 
the periodicity does not alter the function space). On the left hand side (i.e. for xi < 0, 
denoted with the superscript “-”), the Bloch eigenfunctions are harmonic waves and 
the Bloch expansion coincides with a Fourier expansion. We collect properties on the 
left half-domain in Remark 13.61 


Lemma 2.5 (Bloch expansion). Let K E be the number of cells in each direction, 
let u E L^{KYyC) be a function on the sguare {0,Ke) x {0,Ke). Expanding u in a 
pre-Bloch expansion and then expanding each in eigenfunctions we obtain, with 
coefficients E C, 

00 

“W= E E“7.R™(^) 

J&Qk^Qk m=0 

and similarly, for an expansion corresponding to constant coefficients a'^ = 1, 
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To shorten notation, we will use the multi-index A = {j,'rn) in the index-set Ik '■ = 
^ Qk X Qk^ ^ £ No} C I := Z X No- Abbreviating additionally 


U^(x) := 


we may write the formulas of Lemma |2.5| as 

.,±iTt± 


X=(j,m)£lK 


X) 


( 2 . 8 ) 


(2,9) 


The expansion holds for the basis functions with coefficients and for the basis 
functions with coefficients a^. Moreover, due to L^-orthonormality of the functions 
with h = eK and KY^ = (0, h) x (0, h). 


\u\ 


E 




±|2 


( 2 . 10 ) 


3 Outgoing wave condition 

3.1 Poynting numbers and projections 


Let A = (j, m) G I 
Denoting by ei = (1,0) G 


be an index and let 


be the corresponding Bloch function. 
^ the hrst unit vector, we connect to A G / the real numbers 


P/ := Im-/ U^{x)ei ■ [a%x)VU^{x)] dx, 


Ye 


P^ :=Im-/ U^{x) ei ■VU^{x) dx . 

JYe 


(3.1) 


The number P^ is related to the Poynting vector of the Bloch eigenfunction . It 
corresponds (up to a positive constant) to the horizontal group velocity of this eigen¬ 
function and measures its energy flux in horizontal direction: In the case P^ > 0, the 
energy of the wave is travelling to the right, in the case < 0, the energy of the wave 
is travelling to the left. 

Let us point out the relation to Maxwell’s equations: If u denotes the out-of- 
plane magnetic held, i.e.H = (0,0,n), then the electric held is (Pi,P2,0) with 
El = {—iu:e)~^d 2 U and E 2 = {iu:e)~^diu where e is the permittivity of the medium. 
The complex Poynting vector is P = ^E x H, so the real part of its horizontal compo¬ 
nent is Re(ei ■ P) = Re(^Lf 3 P 2 ) = {2u:)~^Be{—iue~^diu) = (2a;)“Rm(-u ocIim), where 
we used that the coefficient a = e~^ is the inverse permittivity. Our expression in (3.1) 


coincides up to the factor 2u with an integral of this expression. Since P represents a 
local energy hux, the physical quantity of a total energy hux is a surface integral over 
P. In fact, for solutions U\ of a Helmholtz equation, the surface integral is independent 


of the position of the surface. Hence our volume integral in (3.1) indeed coincides with 
the physical quantity of a surface integral. 

The index set for X: In our construction, we hx the height h > 0 of the domain 
and the periodicity length e = h/K, the Bloch expansion is performed in this hxed 
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geometry. As a consequence, we consider only indices A = (j,m) G Ir, the frequency 
parameter j must lie in the discrete set Qr x Qr C Z. On the other hand, for arbitrary 
j G Z, we can still consider the functions and U^. They do not depend on A', 
hence also the values are independent of K. 

Definition 3.1 (Index sets and projections). Wc define the set of indices corresponding 
to right-going waves in xi > 0 as 

/J„:={A€/|P+>0}. (3.2) 

The index sets />q, J^q, />q, /<g, J^q are defined accordingly. 

For K ^ N we define the projections II^q as follows: Let u G L‘^{KYy,C) be a 
function with the discrete Bloch expansion 

^\Uf{x). 

XGIk 


Then we set 

n>oM(a;) := 

xeiKOifo 

With this definition, II^q are the projections onto right-going Bloch-waves. The projec¬ 
tions n^Q, n>Q, n<Q, and n^o defined accordingly. 

For k 2 G Qr and I G Nq? the “vertical” projection and the “eigenvalue” 

projection 11^''’ are defined by 

X&{{j,m)&Iji \ j2=k2} 

fir^u{x):= afU^{x). 

X&{(j,m)elK \ m=l} 


Note that the projections of the discrete Bloch expansion indeed coincide with 

of the corresponding vertical pre-Bloch expansion of Dehnition 


2.2 


the projection 11^^ 

The vertical projection is independent of K in the sense that a periodically extended u 


with a larger value of K has the same projection, compare Remark 2.3 


3.2 Bloch expansion at infinity and outgoing wave condition 

We can now formulate the outgoing wave condition for a solution u of the Helmholtz 


equation (1.1). The loose description of our outgoing wave condition (on the right) is: 


The Bloch expansion of u does not contain Bloch waves that transport energy to the 
left. 

For a rigorous dehnition we must deal with the problem that u is not necessarily 
periodic in a:i-direction. Our solution to this problem is to consider u on large domains 
(which reduces the effects of non-periodicity) and to employ a truncation procedure. 
Furthermore, we want to formulate a condition that characterizes u for large values of 
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xi- For these two reasons, we consider u{xi, X 2 ) for xi G {Re, 2Re) with a large natural 
number R» K. 

In order to construct a function on a large domain, we consider the periodic extension 
of u in vertical direction and restrict afterwards the function to a large rectangle. For 
convenience of notation, we restrict our analysis to squares. 

Definition 3.2 (Bloch expansion far away from the interface). Let u G LlM X 
(0,h.);C) he a function on the infinite strip with height h = eK. Let R G NK be a 
multiple of K. We define u : ^ C as the h-periodic extension of u in X 2 -direction. 

We furthermore define functions : RY^ C by 

uj^{xi, X 2 ) := u{Re + xi, X 2 ), (3.3) 

u'i^{xi,X 2 ) := u{-2Re + Xi,X 2 ). (3.4) 

We use the discrete Bloch expansions of the functions G L‘^{RY^; C), 

• (3-5) 

xeiR 

The coefficients {a^pfjx^i encode the behavior of u for large values of |xi|. 

We are now in the position to define the outgoing wave condition for a solution u 
to the Helmholtz equation, using the short notation f := f for averages 

of functions. 


Definition 3.3 (Outgoing wave condition). For K G N, h = Ke, and R G NK, we con¬ 
sider u G X (0, h)]C). We say that u satisfies the outgoing wave condition on the 

right if the following two conditions are satisfied: \u\^ is bounded, independently 

of L >0, and 

(3.6) 




—0 as R ^ 00 . 


RY^ 


Accordingly, we say that u satisfies the outgoing wave condition on the left, if 
fo Il-1 hounded, independently of L <0, and if 


I 12 

(n^Q(M^) —0 as i? —)■ 00 . 


(3.7) 


RY^ 


Let us repeat the idea of condition (|3.6): The function u is considered at the far 


right by construcing u'^ as in Definition 3.2 This function is projected onto the space 


of left-going waves. We demand that the L^-averages of the resulting functions n^Q(M)^) 
vanish in the limit i? —)■ cx). 

With the expansion 


we can write condition (3.6) equivalently as: 
^ 0 as i? —)• 00 . 


E 


\a 


+ 

X,R\ 


(3.8) 


Ae/fli-iz+o 


Our aim is to show that this definition of an outgoing wave condition implies uniqueness 
properties for the scattering problem. 

We note that the uniform L^-bounds for large values of |L| imply, for solutions u 
of the Helmholtz equation, also uniform bounds for gradients, see Lemma A.3 in the 
appendix. 
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3.3 Truncations and m > 1-projections 

In the outgoing wave condition, we study the limit |a;i| ~ i? —)■ cxd and the functions 
on large squares Wr := RY^ = {0,Re)^ with \Wr\ = {eRY- As a measure for typical 
values of a function v we use L^-averages on Wr and the corresponding scalar product, 


{v,w)R:=f v-w 


(3.9) 


Wr 


In the following we denote by Cq = £q = —V ■ (YV) the elliptic operator of (2.7). 
As above, we denote cubes by Wr = RY^ and, by slight abuse of notation, we write 
Wr^i := £(1, R — 1)^ for a smaller cube that has the point £(1,1) as its bottom left 
corner and e{R — 1, i? — 1) as its top right corner. We use a family of smooth cut-off 
functions p '.= rjR with the properties 


G C°°(1Vr; M), riR — lonWR-i, ||V?7r||oo < C'oj IIV^t^rHoo < Co (3.10) 


for some i?-independent constant Cq (e > 0 is hxed), and with compact support in 
(0, Re) X (0, 7?e)tt, where (0, i?e)[t indicates the interval with identihed end points. The 
latter requirement admits sequences rj with compact support in (0, i?e) x (0,i?e), but 
also sequences of vertically periodic functions p, in particular functions r] = r]{xi). In 
the subsequent proofs we do not indicate the /^-dependence of tjr and write only rj. 
We furthermore omit the superscipts ±, the eigenvalue corresponding to A = (j, m) is 
denoted by fix = fimij)- Constants are allowed to depend on e > 0. 


Lemma 3.4 (The effect of truncations). For R E N let rj = tjr be a family of cut-off 
functions satisfying (3.10). Let vr and wr be sequences of functions in L‘^{Wr; C) with 
vr G H^(Wr-,C). We assume that certain averages over boundary strips are bounded: 


If 

R Jwr\Wr_i 


\vr\^ + \Vvr\^<Co, 


1 

R 


'Wr\Wr.i 


\Wr\‘^ < Cq , 


(3.11) 


with Co independent of R. Then, with a constant C that is independent of R: 
1. Application of Co to a truncated function: 


IR-oivR)r] - Co{vRrf)f = 


Wr 


Wr 


Ro{vr)v - X] tYorT], Ux)r Ux 
>^&Ir 


C 

< — . 
- R 


2. If 11 is one of the projections of Definition \3.1\ then 

\U{wr)-I l{wRrf)f < -f \wR-WRrff <^. 


Wr 


(3.12) 


(3.13) 


Wr 


Proof. In the following, the letter C denotes different constants, possibly varying from 
one line to the next, but always independent of R. To prove (3.12), we expand the 
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L^-function Co{vR)ri in Bloch-waves. The following calculation uses several times inte¬ 
gration by parts; due to the r^-factor, no boundary integrals occur. In the hrst equation 
we use that the coefficient a\ in the expansion of C^iy^r] is obtained by taking the 
scalar product with U\ (orthonormality of the U\). 


i^o('>^R)v,Ux)iiUx= ^ {vr,Co{vUx))rUx 






^ {{vrT], CqUx)r + (viiCoiv), Ux)r - 2 {vr a^Vr], VUx)r) Ux 


\eiR 


(hA {vrV, Ux)r - {vrCo^t]), Ux)r + 2 {Wvr ■ a^Vf], Ux)r) Ux 

^&Ir 

/iA {vrT], Ux)r Ux I - VRCo{ri) 2a^VvR ■ Vtj , 

\\eiR J 

where in the third equality we exploited £oh^A = fJ^xUx and fix £ The contribution 
of the last two terms can be estimated by 


/ 

Jwr 


\vRCo{r])f + \2a^VvR ■ Vt]\^ 


— \\^oiv)\\ l°°{Wr)T ^{supp(Vrj)} + || 2 a '^v\\l°°{Wr)T l{supp(V7?)} 


Wn 


Wr 


< 


C 


'Wr\Wr-i 




'Wr\Wr-i 


I I < §. 


In the second inequality we exploited supp(V77) C (Wr \ ITr_i), in the last inequality 
we used the uniform bounds (3.11). This proves the inquality in (3.12). 

Regarding the equality in (3.12) we have to verify that the formal equality Cow = 
PA('iy, f^A)/? f^A holds for functions w G H‘^{Wr) with vanishing boundary data. We 
hnd this from 

{C^w, Ux) = {w, CoUx) = fix{w, Ux) , (3.14) 

where we used in the marked equality that boundary terms vanish. 

Inequality (3.13) is a direct consequence of linearity and norm-boundedness of the 
projections: 



- U{wRr])f 


/ \niwR{l-v))f 

Jwr 


<-[ \wR{l-7])f < 

Jwr 


-J 

Jwr\wr_-, 


\wr? < 


This concludes the proof. □ 

A warning concerning the non-periodicity of truncated solutions. Let u be 

a vertically periodic solution of the Helmholtz equation on M x (0, h) and let u\ be 
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defined as in Definition 3.2 Then u~^ is a solution of the Helmholtz equation on the 
open square Wr = RY^. But is not a periodic solution on the square (since it is not 
periodic in horizontal direction). 

This fact implies that certain formal calculations are wrong: Let u~^ have the Bloch 
expansion for some coefficients a\ (every L^-function posesses such 

an expansion). Then, in general, the following identity fails to hold: 


Co - E «A/:o(f/A+)- 

Ae/fl Ae/fl 

To see this, let us assume that the relation (?) holds. Then the Helmholtz equation 

provides = CqU^ = CoJ^xein^^^x “ 'Zxein^^^o{U+) = 

J2x&iji ^>^hx^x- Uniqueness of the Bloch expansion implies ax{uj‘^ ~ hx) = 0 for every 
A G Ir, hence oa = 0 for every A with fr'l ^ ltY. We conclude that the Bloch expansion 
of contains only contributions from those basis functions with = uj"^. 

This is a contradiction for general solutions u: Let a; be a frequency such that 
is not contained in the discrete set of Let furthermore m be a non-vanishing 

Bloch wave for the frequency u. Then the expansion of u does not use only oa = 0, a 
contradiction. 


Lemma 3.5 (Contributions from energy levels m > 1). Let uj satisfy the smallness 


condition (1.7) of Assumption 1.1. Let u G x (0, h);C) he a vertically periodic 

solution of the Helmholtz equation Cqu = u'^u satisfying the uniform L'^-bounds of 


Definition 3.3 Let rj = tjr be a family of cut-off functions as in (3.10). Then, with a 


constant C that is independent of R: 

/ |n;S(«pr<^ 

Jwr -K 


and 


H 


ev,± 
m>l 


Wr 


( ± 


(3.15) 


Proof. We perform the proof for the superscript “-f”. Relation (3.13) applied to uf 
provides 


Wr 


|nSK) - nSK’))r ^ ■ 


Indeed, by the uniform L^-bounds of u, the condition U ^ < Co with Cq 

independent of R is satished. The above inequality implies that it is sufficient to show 
only one of the two relations in ( 3.15[ ), we show the second. 

We now exploit Assumption 0 Due to ( |1.7[ ), there exists 5 > 0 such that \uj'^ — 
/^Ap > d for all A = {j,m) with m > 1. We therefore hnd 

d-f K>i{'^RV)f = d - h\){uRV,Ux)R\^ 


^ = {j,rn)Gl^ 
i>l 


>1 


< I i^^^R d, Ux)^ - {t\Ur V, Ux)ji\^ = I (^o«) V, Ux)ji - {t\Ur V, Ux) 

X&Ir >^&1r 


Rl 












20 


A. Lamacz, B. Schweizer 


Wr 


^oiuR)v - 


\eiR 


c 

< — 
- R 


In the second line we used that Co{u'^) = ui^W^ holds pointwise almost everywhere in 


Wr. In the last inequality we used (3.12), exploiting the uniform hThpounds provided 
by Lemma A.3 This concludes the proof. □ 


3.4 Other radiation conditions and the sesquilinear form 6^ 

In this section we discuss how our radiation condition simplihes in the case of a homo¬ 
geneous medium. We furthermore show that the radiation condition suggested by Fliss 
and Joly in [16] is formally stronger than our condition: Every solution that satishes 
the condition of [T6| satishes also our condition. Finally, we introduce the sesquilinear 
form 6^, which plays a major role in our proofs. The form 6^ can also be used to 
introduce an even weaker form of the outgoing wave condition. 


In free space, Bloch expansions are Fourier expansions 

Let us study the outgoing wave condition in a homogeneous medium. This is our 
situation for xi < 0, indicated with the superscript In a homogeneous medium, 
the Bloch waves are harmonic waves. This allows to give explicit formulas for some 
quantities and, in particular, for the outgoing wave condition. 


Remark 3.6 (Basis functions and Poynting numbers in a homogeneous medium). The 
functions are e-periodic eigenfunctions of the operator Cf = —(V + 2T\'ii/e) ■ 

(V -|- 2 ^]/e). As such, for fixed j G Z, they are harmonic waves, I ^ ^o} = 

^^ 2 mk-x/e | ^ ^ Morc precisely, for every j E Z and every m E N, there exists a 

wave-number k = k{j,m) E ZZ such that 

= iA-tfR. ( 3 . 16 ) 

Accordingly, for A = {j,m), we have Uf{x) = ^ and the Poynting number is 

r - 27r 

= Im 4 (x)ei ■ VUy^ (a:) dx = — {ki + ji) . 

JYe ^ 

In particular, for the first energy level, A = (j, 0), we find 


k = k{j, 0) G argmin \k j\ , 

k&1‘2 


and thus, for j = (ji,j 2 ) with ji 4 |.- 


£ 


P, 


(jfi) 


ji + argmin \ki + ji 

fciSZ 


ji /orjiG[0,1) 

ji-1 forjiEil,!]' 


(3.17) 


The wave f7(jo) right-going in the sense of Definition 3.1 if and only if ji E [0, | 
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Remark 3.7 (Outgoing wave condition in a homogeneous medium). Let u]^ be as in 
Definition's^ Remark 3.6 implies that, for every index X = {j, m) G Ir, there exists an 
index k = k{\) G with Ufi{x) = = ^ 2 -KiR{k+j).x/R£^ Using the shorthand 

notation /(A) := R{k{X) + j) G the Bloch expansion of ufi can be rewritten as a 
Fourier expansion, 


Ur{X) = 




^\,R^\ 


X = 


Ae/fl 


“a,_r ^ 


27 ril(A)-x/Iie 


By Remark 3.6, the basis function (x) is right-going in the sense of Definition 3.1 if 


and only if IfiX) > 0. This simplifies the radiation condition: The function u with the 


truncations on the left satisfies the outgoing wave condition (3.7) if and only if 


as R ^ oo . (3.18) 

n(A)>o 


We emphasize that, in order to evaluate the radiation condition, the map I : Ir ^ Tfi 
need not be evaluated (we know that it is bijective). By expanding ufi in a classical 
Fourier series, ufi{x) = '^ 1^12 with some coefficients the outgoing 

wave condition (3.18) is equivalent to —)■ 0 as i? —)■ cxd. 


U>o 


Comparison to the outgoing wave condition of Fliss and Joly |[T5j 


We claim that the outgoing wave condition of Fliss and Joly is formally stronger than 
our condition. More precisely: Every solution u that satishes the radiation condition 
(1.5) of [16] satishes also our radiation condition (3.6) (we restrict the considerations 
here to the right side Xi —>■ + 00 ). Indeed, let u be as in (1.5), i.e. a hnite sum of 
right-going Bloch-waves plus an exponentially decaying remainder w^{x). In order to 


check (3.6) for u, it suffices to verify, for each of the hnitely many terms, the smallness 


of its n^Q-projection. The smallness of the projection of w~^{x) is clear because of 
the exponential decay of w~^{x) and the boundedness of the projection operator. The 
smallness of the projection for each of the Bloch waves is shown in the subsequent 
lemma. 


Since the expansion (1.5) contains only vertically periodic waves with frequency uj 


that are outgoing, we restrict our analysis to such Bloch-waves U^. 


Lemma 3.8 (Right-going Bloch waves satisfy the radiation condition). Let iF G N 
denote the number of periodicity cells in vertical direction and let Ufi be a Bloch wave 
with X = (j,m) G Z X No; J 2 & Qk = {0 ,^,...,^}, /i+(j) = u^, and Pfi > 0. 
We impose that the freguency u satisfies Assumption 1.1 (hence m = 0). Then, 
NK 3 R ^ 00 : 


as 


n 




-3 0 , 


RYe 


(3.19) 
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Proof. Step 1: Equivalence of (3.19) and (3.20). By L^-orthogonality of the projections 
we can calculate 


1 = 


RYe 


\iK)P = f HMK)i)\" + f |nj„(W)Sr. 


RY^ 


RYe 


which yields the equivalence of (3.19) and ( |3.20 ). 

Step 2: Proof of (3.19). Let A = (j, m) = (j,0) be as in the lemma. Arguing 


exactly as in Lemma 3.5 we conclude that contributions from energy levels m > 1 are 
negligible, 


RYe 


Consequently, for the weighted L^-norm of n^Q((f/^)^) we calculate 

iny((f/Apar= E i(W)tv>>Rr= 5^ i(W)t.(/x)Rr+o(i) 


RYe 


Ae/flO/ 


<0 




E \((Unt,U>.}P + o(l) 


■m^0,J2^32 


In the last line we exploited that due to j 2 G Qr all scalar products with j 2 ^ j 2 vanish. 
Next we show that there exists a constant C = C{X) > 0 such that 




|{(f^A*)tVA>R| < 


(3.21) 


for all A = (j, 0) G Jij n J^g with j 2 = j 2 - Indeed, a direct calculation analogous to that 
of Lemma |A.1| yields 

n(\ V] 1 _ p27ri0i-ji)R 

with (^(A, A) := . In particular |C'(A,A)| < 1 and therefore 




X _ g27rj(ji-ji)i? 


X _ g27ri(Xi-Xi) 


< 




(3.22) 


We now exploit P£ > 0. The eigenvalue /UqIJ) is simple, hence the Poynting number 
P£ = P^Qj is continuous in the wave number j. One thus hnds a positive constant 
5 = (5(A) > 0 such that 


1 — e 




> 6 


for all A = (j,0) G Jr ft J^g with j 2 = j 2 - Together with (3.22) this yields the claim 


(3.21) with C = 
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With estimate (3.21) at hand we conclude 




RYe 


n<o((t^;^)R)| |A = (j,0) G J/jnJ+o with ja = j 2 | + o(l) 


< — + o(l) —)■ 0 as i? —)■ cxD , 
It 


which was the claim. 


□ 


The sesquilinear form b 


Let r] = rjji he a family of cut-off functions as in (3.10). In view of Relation (3.13) of 


Lemma 3.4, the outgoing wave conditions (3.6) and (3.7) are equivalent to outgoing 
wave conditions for the truncated functions := rj. More precisely, they are 
equivalent to the conditions 


Wr 


^0 and f | (n>o(Mfl,,,)| ^0 as R ^ oo . (3.23) 


Wr 


In the proof of our uniqueness result we will use (3.23) instead of the original conditions 


(3.6) and (3.7). In fact, even a weaker form of the conditions is sufficient and we discuss 


this relaxation in the following. 

Corresponding to the energy flux dehnition in (3.1), we associate to a function 
w G C) on Wr = RY^ the Poynting number 


B+[w) 


:= Im+ w{x)ei ■ [a^(a;)Vi/;(a;)] dx. 

JWr 


(3.24) 


The quadratic expression is dehned analogously, with Y[x) replaced by 1. 


Definition 3.9 (Weaker form of the outgoing wave condition). For iC,i? G N with 


R G KN we consider u G x (0,£iL);C) and as in ( |3.23| ). We say that u 

satisfies the energetic outgoing wave condition on the right, if 

(n<onmLo('W^ ,j)) —)■ 0 as R ^ oo . (3.25) 

Accordingly, we say that u satisfies the energetic outgoing wave condition on the left, if 

(n;on^ro(M;^,,)) as r^oo. (3.26) 


In two respects, the condition (|3.25|) is similar to the condition (|3.23|): the function 

contri- 


u is considered at the far right since only is used. In view of Lemma 3.5 


butions from energy levels m > 1 can be neglected and we consider only n^^Q(M(^^). 


IS 


Furthermore, this function is projected to left-going waves, i.e. only II^qLI^ ko('“R, 
studied. The main difference between the two conditions is that, instead of looking 
at the weighted L^-norm, one demands in (3.25) a decay property for the energy-flux 


quantity Bfi. At the end of this section, we will see that condition (3.23) (together with 


the uniform L^-bounds and the solution property) implies (3.25). 
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The definition of B'^ in ( |3.24[ ) snggests to introduce additionally the (nonsymmetric) 
sesquilinear forms 6^ : L‘^{Wr]£.) x H^{yVR]<C) —)■ C, 

b'^{u,v) := j- u(x) ei ■ [a^(a;) Vn(a:)] dx , 

(3.27) 

b~^[u,v)-.=-j- u{x) Cl ■Vv{x) dx . 

Jwr 

The definition is tailored to calculate energy fluxes. The energy flux of the left-going 
contributions of ^ (in the right half-plane) is quantified by 

= Im / n+oM+^(a;)ei ■ [a^(a;)V(n+oM+^)(a;)] dx . 

JWr 

The connection to is given by 

Pt = B%{Ut) = Im i,± (C/f, Ut) . (3,28) 

Let us collect some properties of the sesquilinear forms 6^. 

Lemma 3.10 (Properties of the sesquilinear form 6^). For R gN, the following holds: 

1. Orthogonality property of 6^. Let X,X^Ir be such that X = {j,rn), X = {j,m) 
with j ^ j. Then U^,U^ of (2.8) satisfy 


bUulU^) = 0. 


(3.29) 


2. Convergence property of 6^. Let sequences ur G L‘^{Wr-, C) and vr G H^{Wr; C) 
be such that 


\ur\‘^ + |Vn_R|^ < Co 


(3.30) 


Wr 


with Cq independent of R. Let either \ur\‘^ —)■ 0 or iVuijp —)• 0 as i? —)• cx). 


Then there holds 


Wr 


-*■ 0, 


(3.31) 


Proof. 1. We prove (3.29) for U^, JJX, the argument for 17^ , U- is analogous. We have 


to show that 


A ’ 




Wr 


dx = 0. 


By definition of and IfX there holds 

U+{x) = 

VUfix) = fvTt(a;) + (i2TT]/e) Tt(x) 

A A A 


^i2nj.x/e 
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with ^-periodic functions \I/^, and Vd'h. Due to the ^-periodicity of and since 
J:J ^ Qr satisfy j j, we can apply Lemma A.l of the appendix, which yields the 
claim. 

2. We show the claim for 6^, the argument for b~^ is analogous. The Cauchy-Schwarz 
inequality allows to calculate 


\bUuR,VR) \ = 


UR{x)ei ■ [a^(a;)VnK(a;)] dx 


Wr 


< lla" 


\ur[ 


Wr 


1/2 


ivn«r 


Wr 


1/2 


—)■ 0 as i? —)■ cx), 


which concludes the proof. 


□ 


Lemma 3.10 shows that the outgoing wave condition (3.23) together with the 
bounds of Dehnition 3.3 imply (3.25). Indeed, by (3.23), there holds 


Wr 


|pr+ rr®'"’’''!?; 


)r< 


Wr 




) I —0 as i? —)■ CX) . 


Moreover, satishes | V <C with C independent 

of R due to Lemma A.3 and Lemma |A.4| of the appendix. Lemma 3.10 provides 

Bl (n+oC^l+(n+^)) = Im6+ (n+oC^+g(n+^), (n+^)) ^ 0 as ^ oo , 


and hence (3.25). 

Remark 3.11 (On the sesquilinear form 5^). Another choice of a bilinear form is 

b~l^{u,v) := ]-+ |-u(a:)ei • [a®(a;) Vn(a;)] —v{x)ei ■ [a®(a;)V'u(a:)]| dx . (3.32) 

^Jwr ^ 

With this choice, the energy flux can be calculated as bef ore, since Im5^(M,n) = 
Im6^(-u,n) holds for every u. The properties of Lemma 3.10 remain also unchanged, 
the only additional reguirement would be an -bound also for ur in (3.30). 

The advantage ofbf^{u,v) is that more orthogonality can he expected for b'^ than for 
b^. Essentially, the bilinear form q of (27) in fWf coincides with hf^ (up to a factor 
2 , our coefficient a®, and the fact that we use a volume integral for the averaging). 
In Theorem 3 of an orthogonality property is shown for q, which resembles our 
orthogonality relation (3.29), stating that orthogonality holds also for X = {j,m) and 
A = {j,m) with m fh. Unfortunately, such an orthogonality is only true for basis 
functions corresponding to the same frequency uj, while our analysis of an interface 
would require orthogonality independent of the frequency. 


4 Bloch measures and uniqueness properties 

Our aim is to show uniqueness properties of the transmision Problem |1.2| with incoming 
wave Uinc and outgoing wave conditions. Following the standard procedure of unique¬ 
ness proofs, we consider two solutions u and u of the problem. Due to linearity of the 
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system, the difference v := u — u satisffes again (1.1). Furthermore, it satisfies outgoing 
wave conditions on the left and on the right according to Definition |3.3 , without any 
incoming wave Ui^c- At this point, we have exploited the triangle inequality: Certain 
projections of u and u tend to zero in a weighted L^-norm, hence also the projections 
of V tend to zero. We can not show that v vanishes (indeed, as explained in the intro¬ 
duction, we expect that there exist nontrivial solutions for vanishing U-mc)- But we can 
show that the functions consist, in the limit R —)■ oo, only of vertical waves. The 
right object to study is the Bloch measure associated with v^. 

We recall that the frequency assumption (1.7) implies that, in the limit R —)■ oo, 
the discrete Bloch expansions of contain only modes corresponding to A = (j, m) 


with m = 0, see Lemma 3.5 


4.1 Bloch measures 


In the definition of Bloch measures we use the space of all Radon measures on the unit 
square Z, which we denote as M(Z). It is the dual of the space of continuous functions 
on Z and accordingly equipped with the topology of weak-* convergence. This means, 
in particular, that every bounded sequence in M(Z) has a convergent subsequence. 

Definition 4.1 (Discrete Bloch measure). Let ur G R^IWr^LI) he a sequence of func¬ 
tions with discrete Bloch-expansions 

ur{x) = afU^{x) , 

Ae/fi 

where = a'^iR) depend on R G N. Given these coefficients, for fixed I G Nq, we 
define the l-th discrete Bloch-measure G j\4(Z) by 

^tR-= I“aP (4-1) 

A=(i,06^H 

where 6j denotes the Dirac measure at the frequency j E Z. 

For Ur fixed, is a non-negative Radon measure on Z = [0,1]^. There holds 


OO 


E 

1=0 




(4.2) 


Our aim is to study the limiting behavior i? —)■ cx) of the discrete Bloch measures r 


± 

l,R- 


Definition 4.2 (Bloch measure). For e > Q, K ^ N and h = Ke, let u be a function 
G L^qp(Mx (0, h); C). We consider a sequence NK 3 R ^ oo. 


u 


We extract := u^rj 


according to Definition \3.^ with a sequence of cut-off functions rj = rjR as in ( |3.10[ ) . 
For I 


G N, 


0) 


let he the discrete Bloch measures associated with 


We say that the measure iz^^ G j\4(Z) is a Bloch measure generated by u if there 
holds, along a subsequence R —)■ oo, in the sense of measures (i.e. weak-*). 
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Relation (4.3) is equivalent to the following: for every test-function 0 G C(Z) on 
Z = [0,1]^ there holds 


E ■>(.»■) I 

A=(i,0e/fl 


J z ' J z 


± 

1,00 


as i? —)■ cxD . 


The methods of Section 

tn- 


4.3 


force us to work with cut-off functions, i.e. with Ur^ and 
with Nevertheless, one may also study the discrete Bloch measures u^r associated 
to (without cut-off function). In the limit R —)• oo, the two measures coincide, 

0 

in the sense of measures. In particular, the Bloch measure generated by Ur^ is inde¬ 
pendent of the choice of p. 

4.2 Energy estimates and consequences for the Bloch measure 


Up to this point (with the exception of Lemma 3.5), our considerations have been 


completely abstract in the following sense: Given a function u G x (0, h); C), we 

have constructed restrictions of u to large boxes, projections of these restrictions, and 
hnally discrete and limiting Bloch measures corresponding to u. Except for regularity 
properties, we have not exploited the Helmholtz equation. In this section, we will derive 
relations that express a physical law: energy conservation. This will eventually lead us 
to the uniqueness properties which are expressed with the Bloch measures. 

The subsequent result states that, while left-going waves on the right vanish by the 
outgoing wave condition, right-going waves vanish by energy conservation. 


Proposition 4.3. Let Assumption 1.1 on oj > 0 be satisfied and let v be a solution 


to the scattering problem (1.1), periodic in vertical direction, satisfying outgoing wave 


conditions on the left and on the right according to Definition 3.3. without incoming 
U\nc = 0. For a seguence of cut-off functions rj = rjR as in (3.10) we consider 


wave, i.e. 


m,ri 


■= VrVr = Y^ 


x&Ir^Xr^X’ c.f. Defimtion\3. 


Then, as NK 3 i? —)■ oo. 


E 

A=(i,0) 
AG/^n/ 




and 


<0 


E 

A=0,0) 

AG/jjn/+ 


Kr\" Px^O- 


(4.4) 


Proof. Step 1: Energy flux eguality. For h = sK and R G NK, we consider the special 
cut-off function f}{x) = TIr^x), dehned for x = {xi,X 2 ) as 


0 


if |xi| < eR , 
if eR < |a;i| < 2eR , 
if ki| > ‘2eR. 


We multiply the Helmholtz equation (1.1) with coefficients a = and solution v by the 
test-function ^{x) v{x). An integration over M x (0, h) and integration by parts yields 
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(no boundary terms appear due to periodicity in X2-direction and compact support): 



[a^ I Vn|^ + v dx^v] = uj^ 



d Ini 


Due to the special choice of 'd and a^{x) = 1 for xi < 0, this equation reads 


— Re ph j•2Re 

/ vdx.v-f 


-2Re Jo 


J Re 



va^dx^v= / / lead'd |n|^ — IVn|^} 


On the left-hand side, we recognize the sesquilinear forms 6^ of (3.27). Because of 
periodicity in X 2 -direction, we may write 

^ [^R ~ ^R ~ 



lead'd |n|^ — -d |Vn|^} 


(4.5) 


Since the right hand side is real, taking the imaginary part of (4.5) yields 


(4,6) 


Relation (4.6) is an energy conservation: The energy flux into the domain from the left 
must coincide with the energy flux out of the domain at the right. 

Step 2: Truncations and (m > l)-waves. We start this part of the proof with an 
observation regarding the cut-off functions; we want to have them in the argument of 
the sesquilinear form. Due to Lemma A.3 and the properties of the cut-off functions 
7] = rjR we have 


„± 


± |2 


Wr 


Tr - ^R,v 


+ ivr± - vi.±,r^ < ^, 


(4.7) 


and therefore, by Lemma [3.10[ 

bi -b% 

= b% (4 - 4) + b% {v%^, 4 - J ^0 as R^oo. 

The energy conservation ( |4.6[ ) therefore implies that, as i? —)■ exo, 

b- (n-J - Im 6+ (n+J ^ 0 . (4.8) 

We next decompose the sesquilinear forms 6^ according to the projections of Dehni- 
tion |3.1[ and suppress the superscript “±” in the projection. We exploit sesquilinearity 
of in both arguments and write 

= ImfoJ (nm>l {'>^R,ri} ^''^R,ri) + (^^^=0 {'>^R,ri) ) ^m>l {'^R,ri)) (4-9) 

+ im6+ (n-=o(4J,n-=oKJ) ■ 
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We want to exploit the smallness of m > 1-contributions of Lemma 3^ The regularity 
result of Lemma A.3 together with the properties of the sesquilinear form of Lemma 


3.10 yield that the hrst term on the right hand side of (4.9) vanishes in the limit as 

y Lemma A.4[ which provides that also the 


i? —)■ oo. For the second term we app 
gradient of small; Lemma 


3.10 


implies 


i>s (nr=o Kr ,}. nr>i t’t„)) ^ 0 as ^ (X., 
i.e. also the second term on the right hand side of (4.9) vanishes in the limit. We hnd 
that, as i? —)• oo, 

(n^=o (4,J > n^=o (4,J) + o(l) • (4.10) 

Step 3: Energy flux and outgoing wave conditions. In this step we decompose 

ImtJ (vgj , n™,„ (aj,,)) as follows: 

imt+(n".,(o+,,),n-.„(o+,)) 

= Im(,+ (nj,n“.„ (4,,). «,)) 

+ Im6+ (nj,n».„ («+J, («+J) 

+ im(,+ (nj,n“.„ («+,), («+,)) (4,ii) 

+ (nj,n-.„ («+„), («+„)) 

= irnij (nj,n“.„ (a+J, (a+J) 

+ (n>onm=o (a^.,)), 

where the last equality holds, since for A = {j,m = 0) G and A = {j,m = 0) G />o 
one always has j 7 ^ j and thus the mixed sesquilinear forms vanish due to orthogonality 
in the wave number, cf. Lemma 3.10 Exploiting the outgoing wave condition (3.23) 
on the right or, better, the weaker expression (3.25), we hnd that the hrst term on the 
right hand side of ( |4.11[ ) vanishes in the limit R ^ 00 . Hence 

Imi,+ (n'„L„(a+,),n^% (a+J) 

= Imiit (a+,) , (a+,)) + o(l) rsR^<x>. 

We emphasize that we only used the energetic outgoing wave condition (3.25) in this 
calculation. 

Combining ( 4.10[ ) with ( 4.12[ ) we hnally obtain, as i? —)■ 00 , 

Im 6 (^ (n>on^=o {'^R,r]) i n^gH^^Q (ur ,,)) + o(l) . (4-13) 

Step 4-' Consequences for outgoing waves. We analyze ( |4.13[ ) further, exploiting the 
discrete Bloch expansion of .V ~ ^\,R^X- 

Im 6 j (H^qH^^q , H^qH^^q 

= I] Y1 FR“L''t 


(4.12) 


x=ijfl)&iRnifg A=(i,o)e/fln/Jo 




\a 


A=(j,o)e/fln/+o 


( 177 .( 7 ,+)= 5 ^ 

A=(i,o)erHnr+p 


a 


\,R 


2 p+ 

w • 
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In the last line we again exploited the orthogonality of the sesqnilinear form 6^ in the 
wave nnmber j, see Lemma 3.10, and the relation (3.28) for . We may therefore 
write (4.13) as 


^ |a+^| 2 p+ + o(l). 

A=(i,o)e/Hn/^(, 

On the left, we hnd similarly 

I«a,rP^a“+ o(l). 

A=(i,o)e/iin/<(, 


The energy relation (4.6) together with the sign properties > 0 for A G and 
Px < 0 for A G /<o allows to conclude (|4.4[) □ 


4.3 Proof of Theorem 1.3 and Theorem 1.5 


We study solutions to the Helmholtz equation (1.1). In order to prove Theorems 1.3 


an d[T ^ we have to check conditions that are satished by the support of Bloch measures 
of solutions. We recall the notation 

J%, = {j 6 Z = [0,1]= I = 0 for A = (i, 0} 

for the index set corresponding to vertical waves. 


Proof of Theorem |1.3[ In the following, we consider solutions to the transmission 
Problem 11.21 We are interested in a function v which is the difference of two solutions 
or, equivalently, a solution to the problem without incoming wave. 


Proposition 4.4 (Solutions in absence of incoming waves). Let Assumption 1.1 on 


u > 0 be satisfied and let v be a solution to the scattering problem (1.1), periodic 


in vertical direction, satisfying outgoing wave conditions on the left and on the right 
according to Definition 3.3, without incoming wave. Let with I G Nq, be Bloch 


measures that are generated by v. Then 


Koo = 0 


for I > 1, 


SUPPll'too) C jfo.O ■ 


(4.14) 

(4.15) 


Proof. We only show the statement for the limiting Bloch measures the argument 
for is analogous. Let rP\ the expansion of the truncated 

solution. Then the corresponding discrete Bloch measures are given by 


A=(y0e/fl 




The case I >1: From (3.15) we know that 


X=(i.m)( 


Tn>l 


, + ,2 C 

l«A,i?l 
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and hence 






\a 


\,R\ 


0 


as -R —)■ oo . 


This shows = 0 for every Z > 1. 

The case I = O; We have to show snpp(z/,fgo) C J=oo- To this end, we consider an 
arbitrary test fnnction 0 G C{Z) with 


snpp(0) C {j e Z\X = {j, 0) G J+o U J+o) • 


The ontgoing wave condition (3.6) and Proposition 14.3| yield, in the limit R —)• oo, 


E 


a 


+ |2 
\,r\ 


—)■ 0 and 


E 




x=(j,o)&iRniXo 


\=ij,o)£iRnn 


In the following we assume that supp(0) fl {j G Z | A = (j, 0) G ^ 0, otherwise the 
proof simplihes. We hnd 

Cl := inf P’a^> 0 . 

36supp(i/)) 

Without loss of generality, we assume 0 > 0 (otherwise we consider absolute values). 
For the limit i? —)■ oo we calculate 

f cj)du^j^= i«A,Rr 0 (j)+ 

J Z 4 ^ 4 _ 

^=(.7,o)e/+|3nifl A=(j,o)ei+pn/fl 

jesupp(0) jesupp(ti>) 

< ll</’l|oo ^ ~ l"A,Rr-fA^ 0- 

A=(i,o)e/Hn/+o A=(i,o)e/Bn/+o 


This shows (4.15) for “+”, since 
argument for is analogous. 


0 with support outside J=oo arbitrary. 


The 

□ 


We next prove that, far away from the interface, solutions to the transmission prob¬ 
lem contain only waves that satisfy the frequency condition and the vertical periodicity. 


on 


Proposition 4.5 (Bloch measures and frequency condition). Let Assumption 1.1 
oj > t) he satisfied and let v be a solution to the scattering problem (1.1), periodic 
in vertical direction, satisfying outgoing wave conditions on the left and on the right 
according to Definition \3.^^ with incoming wave U^c or without incoming wave. Let 
^Ooo ^he Bloch measure to I = 0 that is generated by v. Then 


supp(zaJoo) C [j e z\p^{j) =u?, ja e Z/R:} . (4.16) 

Proof. Remark 

{j G Z\j 2 G TijK}. This implies that also the limit measure are supported on 
this set. 


2.3 


implies that the discrete Bloch measures are supported on 
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In order to show (4.16), it remains to check the frequency condition; we proceed as 
in the last proof. Let 0 : Z —)■ M be continuous and bounded with supp(0) fl {j|/io(j) = 
= 0. Arguing with decompositions of the domain of integration, we can consider 
separately a test-function 0 > 0 with the property 0(j) > 0 ^ ho(j) > and a 
test-function (^ > 0 with the property (f){j) > 0 ^ p-o(j) < The arguments are 
analogous and we consider here only 0 as above. 

By continuity of 0 we hnd some 5 > 0 such that — > 6 for every j G supp(0). 

we 


Our aim is to show that — 0- 

have, as i? —)■ oo. 


By dehnition of the Bloch measure z/, 


0,oo 




a 


X,R\ 


V(J) 


A=(i,o)e/fl 




(4.17) 


A=(i,o)e/fl 


The result = 0 is shown once we prove that the right hand side of (4.17) 

vanishes in the limit i? —)■ cx). In order to show this fact, we recall that the coefficients 
dxR are obtained from a Bloch-expansion of the solution at the far right, i.e. cxxji = 
R- We calculate 


(/io(j) -a;^)|a^,RlV(j) 

'^=(iiO)e/_R 

= 0(j) ajij ^^o{3)Ut)R - U+)r] 

'^=(i)0)er_R 

= E ^ [K„, d:oU+)R - u+)r] 

'^=(i)0)er_R 

<^l0lloo f Y i Y 

\x=U,0)&Ir / \A=(i,0)e/fl 

In this calculation we used the following: (1) formula for (2) the eigenvalue prop¬ 
erty of ?7 a with eigenvalue pa = hm(j), (3) integration by parts without boundary terms 
due to the cut-off function tj, (4) Cauchy-Schwarz inequality. Using orthonormality of 
the basis functions f/^ we obtain 



A=(i,o)e/H 

kJwr / \Jwr / 
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Since satisfies uniform L^-bounds and since „ = uj‘^u^ holds up to a small L 


error, the right hand side of (4.17) is small for large i? > 0. This proves = 0 

and hence (4.16) for “+”. The proof for is analogous. □ 


Proof of Theorem 1^. The difference v of two solutions satishes the outgoing wave 
condition without an incident wave. Theorem 1.3 is an immediate consequence of 
Propositions |4.4| and |4.5[ □ 


Proof of Theorem |1.5| . We now provide the proof of Theorem |1.5[ We therefore 


assume that: Assumption | 1 . 1 | on a; > 0 is satished and u is a solution of the scattering 
problem with incoming wave tljnc, which has the wave number k = (^ 1 ,^ 2 )- In par¬ 
ticular, M is a vertically periodic solution of ( 1 . 1 ) such that u and u — f/mc satisfy the 


outgoing wave conditions on the right and on the left. 

Let be Bloch measures that are generated by the solution u. The frequency 


condition (1.7) is satished and we can therefore use Lemma 3.5 As in Proposition 4.4 
case / > 1, we conclude from (3.15) (and the analogous result for that = 0 
holds for every I > 1. Moreover, according to Proposition 4.5 we have that supp(r'^g^) C 

{j G z\fx^{j) = , js e z/a:}. 


Theorem 1.5 is shown once we verify the following property of the Bloch measure 


■^ 0 , 00 - 


supp(i^ito„) C {j € Z|j2 = fe} U JJd J . 


(4.18) 


Proof of (4.18). We consider the projection of u. This function is again 

a solution of the scattering problem. Indeed, by Lemma A.2 


HL 


one 


has 


Ulfu 


X (0, h);C) with periodicity in the a: 2 -variable, and for arbitrary test functions 
(p e C~(]R X (0, h)) there holds 



{nifu) = f [\cp-a^ nif (Vu) = 



nif (Vv^) ■ Vu 



V (n^f (p) ■a^vu = (p u = 




ip n 


vert 

k2 


where we exploited the orthogonality properties of from Lemma A.l and the 

solution property of u. 

As a consequence, the difference v := u — is a solution of the scattering 

problem with vanishing incoming wave (just as the difference of two solutions in the 
proof of Theorem 1.3[ ). The uniqueness statement of Proposition 4.4 implies: Bloch 
measures (for I = 0 ) that are generated by v have their support in vertical waves, i.e. in 
'^=0,0- 

On the other hand, the Bloch measure of is concentrated on waves with 

vertical wave number k 2 , i.e. in {j G Z\j 2 = ^ 2 }- This follows immediately from the 
fact that all coefficients with 72 7 ^ ^2 in the expansion of vanish. 

Since the Bloch measure of u can have its support only in the union of the supports 
corresponding to and u — the claim (4.18) follows. 

Theorem 11.51 is shown. 
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5 Outlook and conclusions 


Remarks on the existence of solutions 


We give some remarks concerning the existence of solutions to the scattering problem. 
In the end, our radiation condition is “the right one” only if, besides uniqueness, an 
existence result can be shown. 

We formulate the following conjecture: Given a non-singular frequency a; > 0, given 
coefficients a = that are equal to 1 in the left half plane and ^-periodic in the right 
half plane, strictly positive and bounded, given hnally an incoming wave Gjnc as in (1.2) 
(possibly with a condition on k), there exists a solution u of the transmission Problem 

Ol 


The idea for an existence proof is the limiting absorption principle (see e.g. uniEa 
|33] for recent contributions): For a positive artihcial damping parameter 5 > 0, we 
consider the equation 


- V ■ ((1 - i5)a(a;)VM^(a:)) = a;V(a;) (5.1) 

for X G G = M X (0, h). Due to the strictly negative imaginary part of the coefficient 
(1 — i6)a{x), this equation admits a unique solution in the Beppo-Levi space 
as can be shown with the Lax-Milgram Lemma. 

To proceed, two properties must be shown. The hrst is: The sequence satishes 
estimates in some function space, uniformly in 5 > 0. Once this is shown, we can 
consider the distributional limit u of the sequence as 6 ^ 0. As a consequence of 
distributional convergence, the limit m is a solution of the Helmholtz equation with 
coefficients a. 

The intricate part of this approach is to show the second property: The limit u 
satishes the outgoing wave condition. We do not see a straightforward argument that 
yields this condition. 


Our outgoing wave condition in a numerical scheme 


The condition of Dehtition |3.3| might seem inadequate for numerical purposes on hrst 
sight. But it is possible to interpret the condition for a discrete realization: Instead of 
demanding that limits i? —)■ oo vanish, we rather demand that, for some large distance 
R, the left hand side of (3.6) vanishes. In the numerical scheme, this amounts to 
imposing that a hnite number of projections to left-going waves vanishes. 


Conclusions 

We have investigated the transmission properties at the boundary of a photonic crystal. 
Our theorems justify the following: An incoming wave generates, inside the photonic 
crystal, only those Bloch waves, for which the eigenvalue coincides with the (squared) 
frequency of the incoming wave. Furthermore, only those Bloch waves can be generated 
that have the same vertical wave number as the incoming wave; this latter statement 
is true up to vertical waves. 
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Our results rely on a new outgoing wave condition in photonic crystals. The new 
radiation condition is based on Bloch expansions. It is accompanied by a (weak) unique¬ 
ness result, which is expressed with Bloch-measures. The uniqueness result is the basis 
for the analysis of the transmission problem. 
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A Orthogonality and regularity properties 

Lemma A.l (Orthogonality with periodic weight). Let f : M —>■ C be e-periodie and 
integrable, let R eN be an integer. 

1. Orthogonality of exponentials. Let j,j G Qr with j ^ j. Then 

psR 

/ = 0 . (A.l) 

Jo 

2. Orthogonality of the vertical pre-Bloch projection. Let u,v E T^q^(M x (0, ei?); C) 
and let k 2 E Qr. Then there holds 

peR _ peR _ 

/ f{y)u{.xi,y)fllfv{xi,y)dy= / f{y)Y\lfu{xi,y)Y[lfv{xi,y)dy. 

•Jo Jo 

(A.2) 


Proof. 1. By dividing the interval (0,£i?) into subintervals of length e, we obtain 


rsR 


R-l 

fyy2nijy/Se-2nfjy/s dy = 

k=0 ■ 


'‘{k+l)s 


fiy)e 




dy 


' ke 


i? 1 

/ f{y + ke) 

k=0 


R-l 


ym{j-j){y+ke)/e 1 g27ri(j-i)fc 


f{y)e 


2Tri{j-j)yl£ 


dy, 


k=0 


where in the last equality we exploited the periodicity of the weight /. By setting 
C{j,j) ■= Jq dy we conclude 


f'sR 


R-l 


f{y)e 


'^TTijyle -2mjyl£ j _ 


dy = c{j,j)Y^ 


k=0 




= C{3,3) 


X _ g27riO-X)K 
X _ Q2ni{j-j) 


= 0 . 


In the last step we used j,j E Qr, which implies R{j — j) E Z and j,j < 1, and 
exploited j j. 

2. Let u, V have vertical pre-Bloch expansions 

u{xi,X2)= ^ v{xi,X2) = ^ ^-^(a;!, 0 : 2 ) 

L&Qr 
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Then the left hand side of (A.2) reads 

peR 


'0 


f{y)u{xi,y) Ulfv{xi,y)dy 


reR 


/ fiy)^n (^ 1 . y) (^ 1 , y) dy . 

32 &Qr 


Since the function •)$fc 2 (a;i, •) is e-periodic, we can apply the orthogonality 

) of Item 1. The sum on the right hand side collapses to j 2 = ^2 and we hnd 


Lemma A.2 (Vertical pre-Bloch projection and gradients). Let A G N, h = eK, and 
^2 ^ Qk- Let u G X (0, h); C) he periodic in the X 2 -variable. Then the function 

G X (0, h); C) is periodic in X 2 and there holds 

V (ng'V) = ng" (Vt.). (A.3) 

Proof. Let u have the pre-Bloch expansion u{xi,X 2 ) = XljaeQif^ 2 ) 

Due to the periodicity of u in the X 2 -variable, each in the above (hnite) sum has 
A^-regularity, and thus 


(A.2 


Vu{xi,X2)= 

j2&QK 

= ^ [W^j^{xi,X2) + 2'Kii2/e^j.,{xi,X2)e2]e^'^'‘^^'^^/'', 
j2&QK 


(A.4) 


where 62 = (0,1) G denotes the second unit vector. Since the expression in the 
squared brackets is ^-periodic, (A.4) is an expansion of Vm; uniqueness of the pre- 
Bloch expansion implies 

n^f (Vm) {xuX2 ) = {V<!>k,{xi,X2) + 2mk2/e<l>k,{xi,X2)e2)e^-^’^^^^/^ 

= V {<!>k,{x^,X 2 ) = V {Ulfu) {xuX 2 ), 


which proves (A.3). 


□ 


Lemma A.3 (Caccioppoli estimate). Let u G x (0,h)) he a vertically periodic 

solution of the Helmholtz equation Cqu = u'^u. Let u satisfy the uniform L'^-hounds of 
Definition\3.3. Then there holds 


1 

R 


\u 


'Wr\Wr-i 


±\2 

r\ 


+ |Vm^|^ < C and 


\u 


±|2 

r\ 


+ iv4r' < c 


(A.5) 


Wr 


with C independent of R. 

Proof. The proof is, up to translations and a summation, analogous to the proof of the 
standard Caccioppoli estimate: On a rectangle (L — 1,L -|- 2) x (0, h) we use a cut-off 
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function 0 with compact support that depends only on xi and which is identical 1 on 
(L,L + 1) X (0, h). Testing the equation with provides 


r-L+2 I>h 


i‘L+2 ph 


2 | | 2/)2 
u \u\ 0 = 


Cqu{9‘^u) = 


nL+2 ph 


lL -1 Jo 


IL -1 Jo 


'L -1 Jo 


{a^\Vu\^e^ + 2a%Vue) ■ (V0h)} 


The Cauchy-Schwarz inquality is used to treat the last term, the hrst factor is absorbed 
with Young’s inquality in the gradient term, the other consists (up to bounded factors) 
only of the L^-norm of u. We conclude that a bound for the L^-norm on {L — 1,L + 
2) X (0, h) implies a bound for the L^-norm of the gradient on {L,L + 1) x (0, h). A 
summation over many squares yields the result. □ 

Lemma A. 4 (Regularity of eigenvalue projections 11®''). Let be a sequence of 

functions with -regularity and vanishing boundary data, i.e. vr G Hq(Wr;C). We 
assume that 


|^ij|' + |Vu«|V|£o(^i?)|'<^^o 


(A.6) 


Wr 


holds for Cq = —V ■ (a*^ V) with some R-independent constant Cq. 


1. Let n be any of the projections of Definition 3J_. Then there exists an R- 
independent constant C such that 


Wr 


|V (nS<'R)t + |v (n“iy.R)t + |V (n (nS-R))!' < c. (a.t) 


2. If, additionally, —)■ 0 as R ^ oo, then there holds 


IWr 


Wr 


V 1^0 as R^oo. 


(A.8) 


Proof. 1. We omit the superscripts ±. Concerning (A.7) we note that, because of 
= Vr — the estimate for follows directly from the estimate for 

and Assumption (A.6). 

Since = J2x={j o)eiR is a hnite sum of periodic functions, we hnd that 

is periodic in Wr. This allows to calculate, with 0 < a* < inf a^. 


Wr 

( 1 ) 


< 


< 


|v (n^%u«)r < f a^v (n^%u«) ■ V (n-=o^^^) 

'Wr 


( 2 ) 


^0 (n^^gUfi) n^^gUR - -f {Covr) n^^gu/? 


Wr 


Jw 


Wr 


Wr 


|^®J=o(£o^^«)| 

.- . 


1/2 




Wr 




Wr 
1/2 


1/2 


Wr 


<Co. 
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In (1) we exploited the periodicity of to perform integration by parts withont 

boundary terms. In (2), we used the periodicity of vr, which yields £o (n^^gUj?) = 
{CqVr), as shown in (|3.14[). In the last line we exploited the norm-boundedness 


of projections. The claim for 11 (n^^gU/^) is shown analogously, using again periodicity. 
This concludes the proof of Relation ( |A.7 ). 

2. The proof of Relation (A.8) is similar and can be interpreted as an interpola¬ 
tion between function spaces. Once more, we exploit that vr has vanishing (and thus 
periodic) boundary data and that is periodic as a hnite sum (see Item 1.). 

Therefore also the difference If^^^nR = vr — is periodic. Arguing as above we 

obtain, as R —)• oo. 


/ |v(n->i!'R)f < (/ 

Jwr \Jwr 


\^oVr\‘ 


1/2 




Wr 


1/2 


^ 0 . 


(A.9) 


This shows (A.8) and concludes the proof. 


□ 
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